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Abstract 

Concerning the problem of shape and pattern de- 
scription, a compact formula (GiehV formula 1-3 ) has 
recently been proposed that generates a vast diversity 
of natural shapes. However, this formula is a modi- 
fied version of the equation for the circle, so that it is 
expressed in terms of trigonometric functions which 
are inherent to linear phenomena but rarely appear 
in the description of nonlinear phenomena. In conse- 
quence, it generates highly idealized (Platonic) forms 
rather than real- world forms. Since natural shapes 
and patterns generally appear as a result of nonlin- 
ear dynamical processes 4-7 , it would be natural that 
formulas which aim to describe them should be ex- 
pressed in terms of related nonlinear functions. Here, 
two examples of simple mathematical formulas which 
are natural nonlinear modifications (one being a gen- 
eralization) of Gielis' formula are discussed. These 
formulas involve a comparable number of parameters 
and provide non-Platonic representations of a vast di- 
versity of natural shapes and patterns by incorporat- 
ing diverse aspects of asymmetry and seeming disor- 
der which are absent in the original Gielis' formula. It 
is also shown how diverse sequences resembling some 
natural-world pattern evolutions are also generated 
by such nonlinear formulas. — 

Although nonlincarity is an ubiquitous feature of 
real-world phenomena, some aspects of this funda- 
mental property remain surprisingly poorly explored, 



partly because they are rarely incorporated into the 
mathematical tools and models which aim to describe 
even the simplest of such phenomena. In the context 
of the description of natural forms, such as biolog- 
ical shapes, a recent example is provided by Gielis' 
formula 1 : 
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where (r, 9) and (a, b, c, m, ri2, n.3) are polar coordi- 
nates and parameters, respectively. While formula 
(1) generates a certain variety of abstract symmetric 
shapes, which in some cases are reminiscent of nat- 
ural forms, it is clear that its linear origin prevents 
it encompassing most natural forms, these being of- 
ten characterized by diverse types of asymmetry and 
different degree of apparent disorder. Clearly, this 
is because Gielis' formula essentially describes modi- 
fied circles, and hence the use of harmonic functions 
in its formulation. Perhaps the simplest nonlinear 
generalizations of formula (1) would be based on the 
substitution of the trigonometric functions by Jaco- 
bian elliptic functions 8 (JEF's). This choice is quite 
natural since the solutions of the most ubiquitous 
nonlinear integrable oscillators 9 (such as the pendu- 
lum, the Duffing, or the Hclmholtz) are naturally ex- 
pressed in terms of JEF's. This is also the case for 
the most universal nonlinear integrable partial dif- 
ferential equations, 10 ' 11 such as the sine-Gordon, the 



Korteweg-de Vries, the nonlinear Schrodinger, or the 
nonlinear heat equation. 

Here the following two nonlinear (elliptic) modifi- 
cations of Gielis' formula are proposed: 
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where am(«; ra) is the JEF of parameter m, K(m) 
is the complete elliptic integral of the first kind, 
and Up, tp' , A, A') are additional parameters. Note 
that ri{9) reduces to r(9) for m = tp = tp' = 0, 
while m (9) only coincides with r (9) at 9 = for 
A = A' = ip = tp' = 0. An initial shape is generated 
over the 9 interval [0, 2ir]. The dependence of the 
elliptic parameter on the angle coordinate in r H (9) 
is the key feature providing diverse variations of a 
given initial shape. The rate and nature of such 
variations on an initial theme (pattern) can be 
controlled by changing the parameters. Thus, one 
can obtain sequences that mimic transformations of 
biological shapes, including growth processes, and 
evolutions of vortex boundaries in fluids 12 , inter 
al. The top row of Fig. 1 shows an illustrative 
example of this property for the case of a schematic 
(daisy-like) flower whose initial ellipse-like petals 
transform into triangle-like ones, ft is worth men- 
tioning that even such simple representations of 
flowers cannot be generated by formula (1). To 
achieve that goal, r(9) would have to be multiplied 
by an additional harmonic function 1 . The sequence 
was generated using the type 11 formula at the 



parameter values (a, b, c, n\, ri2, n%, ip, <p', X, A ) = 
(2,173,12,-2,2,2,0,tt,1,1) over the 9 intervals 
[0,2tt], [37t,57t], [207r,227r], and [9220tt, 9222vr], 
respectively. Another important feature of real- 
world shapes that is incorporated by the above 
nonlinear formulas is a certain disordered-looking 
impression, such as that in a rose after a spring 
shower. This feature is illustrated in the second 
row of Fig. 1, where different versions of a rose- 
like flower are depicted. The first is generated 
using Gielis' formula at the parameter values 
(a, b, c, n\,ni, n 3 ) = (2, 5, 1.61803, 5, 5, 5), and would 
represent an idealized symmetric rose. This con- 
trasts with the second and third versions, which 
were generated using the type II formula at the 
additional parameters Up,tp',X,X') = (0, it, 1,1) 
and (0,0,0.995,0.995), respectively, and with the 
fourth version which was generated using the 
type I formula at the additional parameter val- 
ues (<p,ip',m) = (0, 7T, 0.999). The third row of 
Fig. 1 shows a sequence of starfish-like patterns 
to illustrate the flexibility of the above nonlin- 
ear formulas in reproducing different aspects of 
asymmetry. Gielis' formula generated the perfectly 
symmetric first version at the parameter values 
(a, b, c, rii, ri2, ns) — (2,5,10,3,3,3). The type II 
formula generated the second and third versions at 
the parameter values (a, b, c, n\, rii, n^, tp, tp', A, A') = 
(2, 5, 10,3,3,3,0,7T, 1,1) and (2, 10, 10, -2, 2, 
2, 0, 7r, 0.91, 1), respectively, and the type 

1 formula generated the fourth version at the 
parameter values (a, b, c, n\, rii, n%, tp, tp' , ra) = 
(2, 10, 10, -2, 2, 2, 0, 7T, 1 - 10~ 14 ) . Formulas 
types I and II are examples of an elliptic class 
that should certainly be useful in graphics and 
modeling where several basic shapes may com- 
bine to yield more complex pictures. Figure 

2 shows a flower obtained by superposition of 
two type I shapes generated at parameter values 
(a, b, c, m, n,2, ns, tp, tp' , m) =(2, 3, 10, 3, 3, 3, 0, tt, 0.95) 
and (2, 17, 10, 3, 3, 3, 0, tt, 0.999). 

In conclusion, since the price of the strikingly broad 
and complex diversity of natural shapes is the inher- 
ent nonlinear character of the processes giving rise 
to them, one should expect that even the most ele- 
mental formulas aimed at describing them should be 
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expressed in terms of nonlinear mathematical func- [12] Friedland, L. & Shagalov, A. G. Resonant forma- 
tions closely related to general (i.e., nonlinear) natu- tion and control of 2D symmetric vortex waves, 
ral processes. The formulas presented here are solely Phys. Rev. Lett. 85, 2941-2944 (2000). 
intended as examples of this line of thinking. The 
challenge now is, as in the case of Gielis' formula, 
to relate them (or similar formulas) to some basic 
natural mechanism of self-organization and pattern 

formation. 0.0.1 Figure Captions 
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Figure 1. Shapes generated using nonlinear formulas 
(2)-(6). 

Figure 2. Flower composed of two shapes gener- 
ated by formula type I. 
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